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I.  INTRODUCTION 


The  Analytic  Sciences  Corporation  (TASC)  recently  performed  a  study  (see 
Reference  1)  for  MICCM  on  the  Implementation  of  a  strapdown  Inertial  updated 
midcourse  guidance  scheme  for  application  to  Short  Range  Air  Defense  type 
missiles.  Part  of  this  effort  Involved  design  of  an  Extended  Kalman  Filter 
(EKF)  to  process  the  radar  measurement  data  received  from  a  "quiet  radar," 
i.e.,  track  on  scan  radar,  as  it  scans  a  target.  The  target  system  dynamic 
model  was  assumed  to  be  linear  in  a  rectangular  coordinate  system  and  the 
target  state  vector  was  defined  as 

x  -  (x,  y,  z,  x,  y,  z,  X,  y,  2) ,  (1 

where  (x,  y,  z)  denote  position,  (x,  y,  z)  denote  velocity,  and  (X,  y,  X) 
denote  acceleration.  Since  the  target  accelerations  cannot  be  directly 
measured,  they  are  modeled  as  first-order  Markov  processes  given  by  the  dif¬ 
ferential  equations 

x  -  -Vjjk  +  w* 

7  “  -*yy  ■*' {Wy  '  '  (2 

■  }  i  , 

z  -  -\zir  +  wz 

where  Vj,  \y,  X2  are  specified  bandwldths  (nominally  equal  and  ^  0)  and 
Wj,  Wy,  wz  are  white  Gaussian  noise  processes  with  zero  means  and  specified 
spectral  densities. 

The  dynamic  system  model  for  the  target  tracking  filter  which  processes 
the  radar  measurement  data  was  given  (see  Reference  1)  as 

x  -  Fx  +  w  (3; 

where  F  is  a  time— invariant  9x9  matrix  and  w  is  the  system  process  noise  vec* 
tor.  Expressing  equation  (3)  in  expanded  form  results  in 


In  order  to  obtain  a  solution  for  x  from  Equation  (3)  an  expression  must  be 
found  for  the  system  state  transition  matrix  defined  by 

®  -  e*At,  (5) 

where  At  is  the  time  interval  (assumed  to  be  constant)  between  measurements 
and  F  is  defined  as  in  Equation  (4). 

There  are  several  methods  (see  Reference  2)  which  may  be  used  to  evaluate 
a  matrix  exponential  of  the  form  eFt.  The  two  more  commonly  used  are: 

(a)  expansion  of  e?t  into  a  power  series  in  t  aaf 

eFt  .  i  +  Ft  +  +  -  ,  (6) 

and  (b)  the  Laplace  transform  method  whereby  eFt  is  found  as 

eFt  -^“l{(si  -  F)-l}.  (7) 

As  the  order  of  the  F  matrix  increases,  the  computations  involved  in  obtaining 
an  exact  evaluation  of  e?t  using  Equation  (7)  become  extremely  burdensome.  As 
a  consequence,  TASC  chose  to  approximate  e'^t  by  using  the  first  two  terms  of 
a  power  series  expansion  in  At,  l.e., 

eFAt  a  j  +  FAt>  (8) 

which  is  very  simple  to  calculate. 

An  exact  expression  will  be  developed,  in  this  report,  for  the  9-state 
target  model  state  transition  matrix  defined  by  Equation  (5)>with  the  matrix  F 
of  Equation  (4).  In  the  process,  another  method  (based  on  the  LaGrange- 
Sy Ives ter  interpolation  polynomial  (see  References  2  and  3)  which  can  be  used 
to  evaluate  eFt,  and  which  enables  one  to  arrive  at  an  exact  expression  for 
e?t  without  requiring  matrix  inversion,  will  be  illustrated. 

II.  MATHEMATICAL  BACKGROUND 

The  method  to  be  used  in  evaluating  e^  will  be  presented  in  this  section, 
and  an  example  will  be  worked  to  illustrate  application  of  this  method  in  com¬ 
parison  with  the  Laplace  transform  method.  First  though,  some  terms  which 
must  be  defined  are  as  follows: 

(a)  Annihilating  Polynomial:  An  annihilating  polynomial  of  a  square 
matrix  F  is  a  scalar  polynomial  f(\)  for  which  f(F)  «  0. 

(b)  Monlc  Polynomial:  A  monic  polynomial  is  a  scalar  polynomial 

f(\)  ■  +  ....  +  cfc  in  which  cQ  ■  1. 

(c)  Minimal  Polynomial:  A  minimal  polynomial  of  a  matrix  F  is  the  monlc 
annihilating  polynomial  of  least  degree. 

(d)  Spectrum:  The  spectrum  of  a  matrix  F  is  the  set  of  characteristic 
values  of  F,  i.e.,  the  set  of  all  scalars  \  for  which  F-\I  is  not  invertible. 


A  function  f(\)  is  said  to  be  defined  on  the  spectrum  of  a  matrix  F  if  the 
values  of  the  function  f(X)  on  the  spectrum  of  F  exist,  i.e.,  if  the  m  numbers 


f(^k).  f’(Mc).  -  .  f(nk"1)(Ak)  (k  “  1.  2, 


8) 


have  meaning,  where  s  is  the  number  of  distinct  characteristic  values  of  F  and 
ntfe  is  the  multiplicity  of  the  k-th  characteristic  value  in  the  minimal  poly¬ 
nomial  (see  Reference  3). 


The  Lagrange-Sylvester  interpolation  polynomial  for  a  function  f(X) 
defined  on  the  spectrum  of  a  matrix  F  is  a  polynomial  r(X)  defined  as 


r(\) 


/  ®k 

[l 

1  \j-l 


akj(^  " 


(9) 


where 


akj 


(j-l)J 


fete)] 


(j  -  1.  2, 

IX  -  Xk 


•.  «k»  k  -  1,  2, 


s),  (10) 


¥k(X)  is  the  minimal  polynomial  of  F  with  the  Xk  term  removed,  and  r(X) 
satisfies  the  conditions  r(Xk)  ■  f(Xk),  r'(Xk)  » f'(Xk),  — ,  r^mk  ^(^k)  “ 

f(mk  15 (Xk),  (k  -  1,  2,  - ,  s).  The  polynomial  r(X)  assumes  the  same  values 

on  the  spectrum  of  F  as  does  f(\)  and  one  has  f(F)  -  r(F). 


If  the  equation  for  the  coefficients  (Equation  (10))  is  expanded  and 
substituted  into  Equation  (9),  Equation  (9)  can  be  expressed  (see  Reference  3) 
in  the  form 


s  mk  .  . 

r(X)  -  l  l  fU  ;(Xk)$kj(\)  (11) 

k-1  j-1 

where  the  (X)  represent  polynomials  in  X,  of  degree  less  than  m,  which  are 
completely  determined  when  ?(X)  is  given  and  are  independent  of  the  choice  of 
f(X).  Equation  (11)  can  be  rewritten  to  obtain  the  fundamental  formula  for 
f(F)  (see  Reference  3): 

f(F)  -  I  I  f<J-1)(Xk)Zkj  (12) 

k-1  j-1 

where  Zkj  -  $kj(F)  (j  -  1,  2,  - ,  n^;  k  -  1,  2,  — ,  s).  As  for  the 

<hkj( X)  terms  of  Equation  (11),  the  Zk<  matrices  are  completely  determined  by  F 
ana  are  Independent  of  the  choice  of  f(X)>  The  Zkj  matrices  are  called  (see 
Reference  3)  the  components  of  the  given  matrix  F. 

There  are  several  methods  by  which  the  Zkj  matrices  can  be  evaluated. 

The  one  which  will  be  Illustrated  here  is  the  method  whereby  certain  simple 
polynomials  for  f(X)  are  successively  substituted  into  the  fundamental  formula 
given  by  Equation  (12).  From  the  linear  equations  so  generated  the 
Zkj  matrices  can  be  determined.  For  comparison,  the  Laplace  transformation 
method  for  finding  e^*  will  also  be  illustrated. 


r,  TV*  '  v'-  '  v'  v' 


V"\  H  ». 


Example  1:  Lagrange-Sylveater  Method 


Suppose  one  wishes  to  find  an  expression  for  eFt  where  F  Is  given  as 

fo  1  ol 


F  - 


2 


-5 


(13) 


For  this  given  matrix  F  one  has  A(X)  -  ¥(X)  ■  (X-l)2  (X-2),  thus,  X^  ■  1, 

\2  ■  2,  »i  -  2  and  m2  ”  1.  Checking  to  see  that  f(X)  -  e^c  is  defined  on  the 
spectrum  of  F  it  is  seen  that,  for  the  root  X^  *  1  with  m^  -  2, 


f(X^)  ”  ec  exists  and 


f' (X^)  ■  ec  exists. 

Also,  for  \2  ■  2  with  m2  ■  1  it  is  seen  that  f(X-2)  *  e^t  exists.  Hence, 
f(X)  ■  e^c  is  defined  on  the  spectrum  of  F. 


Next,  making  use  of  Equation  (12),  f(F)  can  be  written  as 

f(F)  .  eFt  -  f(l)Zu  +  f’(l)Zi2  +  f(2)Z2l  (14) 

and  the  components  of  F  can  be  determined  as  follows.  First,  take  r(X)  to  be 
r(\)  ■  (X-l)2.  For  this  choice  of  r(X),  one  has 

r(F)  -  (F-I)2  -  Z21.  (15) 

If  one  next  takes  r(X)  as  r(X)  »  (X-l)(X-2)  then  one  has 


r(F)  -  (F-I)(F-2I)  -  (-1)  Zl2 
and ,  therefore , 

Z12  -  -  (F-I)(F-2I)  .  (16) 

Finally,  take  r(X)  as  r(X)  ■  1.  In  this  case,  one  has 
r(F)  -  I  -  ZU  +  Z2i 
and  Zn  is  found  to  be 

Z11  "  I  -  Z2i  -  I  "  (F-I)2.  (17) 

If  Equation  (15),  (16)  and  (17)  are  substituted  into  Equation  (1A),  the 
resultant  expression  for  is  found  as 

eFt  ■  e^^Zjj  +  te^^tZi2  +  ®^^tz21 

-  [I-(F-I)2]ec  -  [ ( F-I ) ( F-2I ) ] te c  +  (F-I)2  e2*  • 


(18) 


Performing  the  matrix  multiplication  indicated  in  Equation  (18)  will  result  in 
the  following  expression  for  eFt: 


»2t 


“  o  2  -r 

"2  -3  l" 

"l  -2  l“ 

ePt  - 

-2  5  -2 

et  _ 

2-3  1 

te*  + 

2-4  2 

_-4  8  -3_ 

1 

H 

? 

M 

_ 1 

- 1 

GO 

1 

r**2tet  +  e2t  (2+3t)ec  -  2e2t  (-l-t)et  +  e2t  ' 

(-2-206*  +  2e2t  (5+3t)e*  -  4e2t  (-2-t)e*  +  2e2t  |  .  (19) 

|_(-4-2t)e*  +  4e2t  (8+3t)e*  -  8e2t  (-3-t)e*  +  4e2tJ 

Example  2:  Laplace  Transform  Method 

To  evaluate  e**  using  the  Laplace  transform  method  given  by  Equation  (7), 
the  sequence  of  operations  to  be  performed  would  be  as  follows.  First,  form 
the  matrix  (sI-*F)  as 


(sI-F) 


s 

0 

-2 


-1 

s 

5 


0 

-1 

s-4 


The  determinant  of  (sl-F)  is  A(s)  ■  (s-1)2  (s-2).  The  matrix  inverse, 
(sI-F)~l ,  must  next  be  calculated  and  is  found  to  be 


(sI-F)-l 


(s-1)2(b-2) 


,2  _ 


L  2s 


+  5 


s-4 
s(s-4) 
-5s  +  2 


(20) 


The  expression  for  e**  is  calculated  by  taking  the  inverse  Laplace  transform 
of  each  term  of  (sI-F)“l.  If  Equation  (20)  is  expanded  using  partial  frac¬ 
tions,  the  result  is  given  by 


(sI-F)-l  - 


2 

.  J_ 

3 

2  _ 

2 

/_  1 

.  _1_  + 

1  > 

(s-l)2 

s-2 

(8-1)2 

T 

s-l 

s-2 

V  (8-1)2 

s-l 

s-2^ 

2 

3 

5  _ 

4 

L  i 

_  _2_  + 

2  \ 

(s-l)2 

8-1  s-2 

(8-1)2 

T 

s-l 

s-2 

'  ( 8—1)2 

s-l 

s-2' 

2 

3 

8  _ 

8 

L  1 

_  J_  + 

4  \ 

(s-l)2 

8-1  8-2 

(s-l)2 

s-l 

s-2 

v  (8-1)2 

s-l 

s-2' 

(21) 


and  taking  the  Inverse  Laplace  transform  of  each  term  then  gives  the  final 
expression  for  e**  as 


-Ite*  +  e2t  (3t+2)ec  -  2e2t  (-t-l)e*  +  e2t 
(-2t-2)et  +  2e2t  (3t+5)ec  -  4e2t  (-t-2)e*  +  2e2* 
(-2t-4)ec  +  4e2*  (3t+8)e*  -  8e2*  (-t-3)e*  +  4e2t 


(22) 


,Ft 


J 


Comparing  Equations  (22)  and  (19),  it  is  seen  that  both  approaches  result  in 
the  same  answer;  however,  evaluation  of  eFt  using  the  Laplace  transform  method 
involves  more  of  a  computation  burden  than  does  the  method  using  matrix  com¬ 
ponents,  even  for  a  3x3  matrix. 

Ill.  TARGET  MODEL  STATE  TRANSITION  MATRIX  EVALUATION 
A.  Introduction 

In  this  section  the  matrix  exponential  e^c,  with  F  as  defined  in 
Equation  (4),  will  be  evaluated  to  obtain  an  exact  expression  for  the  target 
state  transition  matrix  for  several  different  combinations  of  X*,  Xy,  \z.  The 
first  step  in  the  process  is  to  evaluate  the  determinant  of  (F-XI)  as 


0 


0 


0  -Xy -X  0 

0  0  -x2-x 


X6[X3  +  (Xx+Xy+Xz)X2  +  (XxXy+XxXz+XyXz)X  +  XxXyXz]. 


(23) 


The  characteristic  values,  therefore,  are  given  by 

Xi  -  X2  -  X3  -  X4  -  X5  -  Xg  ■  0,  X7  -  -Xx,  Xg  -  -Xy,  X9  -  -Xz.  (24) 

The  next  step  is  to  determine  the  minimal  polynomial  of  F.  The  minimal  poly¬ 
nomial  of  a  matrix  will  contain  each  of  the  distinct  characteristic  values  of 
the  matrix  with  a  multiplicity  not  greater  than  each  root  has  in  the  charac¬ 
teristic  polynomial.  The  candidate  choices  for  the  minimal  polynomial  of  F 
are  thus: 


(X+Xx)(X+Xy)(X+Xz)X 

(X+Xx)(X+4)(X+Xz)X2 

!  (25) 

I 

(X+Xx)(X+Xy)(X+Xz)x6 


and  the  actual  minimal 


polynomial  will  depend  on  the  values  chosen  for 


B.  Case  With  Xx,  Xy,  \z  Positive,  Unequal 

Using  each  of  the  expressions  In  Equation  (25)  in  turn,  ¥(F)  is  eval¬ 
uated  to  determine  which  expression  will  result  in  Y(F)  =  0.  Upon  performing 
this  procedure,  the  minimal  polynomial  (for  the  X's  of  this  section)  is  found 
to  be 


¥(X)  -  X2(X+Xx)(X4Xy)(X+Xz). 


(26) 


The  fundamental  formula  for  f(F),  as  given  by  Equation  (12),  will  therefore  be 


f(F)  -  f(X1)Z11  +  f»(Xx)Z12  +  f(X2)Z2l  +  f(X3)Z31  +  f(X4)Z41 


-  f(0)Zu  +  f‘(0)Zi2  +  f(-Xx)Z21  +  f(-Xy)Z31  +  f (-Xz)Z4i .  (27) 

For  f(X)  *  e^c  one  has  f(F)  ■  eFt  so  that,  using  Equation  (27), 

eFt  ”  e^^^Z^j.  +  te^*Z32  +  e^2*Z2i  +  e^3^Z3i  +  e^^tZ4j 

-  (Zu+tZ12)  +  e"XxtZ21  +  e'Vzs!  +  e_KztZ41.  (28) 

The  components  of  F  can  now  be  found  by  choosing  simple  polynomials  for  r(X). 

This  process  is  illustrated  in  the  following  steps  for 

r(X)  -  r(X1)Zn  +  r' (X^)Zj2  +  r(X2)Z2^  +  r(X3)Z3i  +  r(X4)Z41. 

1.  r(X)  -  X(X+Xx)(X+Xy)(X+Xz) 

r'(X)  -  4X3  +  3 (Xx+Xy+Xz ) X2  +  2(XxXy+XxXz+XyXz)X  +  XxXyXz 

F(F+XxI)(F+XyI)(F+XzI)  -  (0)Zn  +XxXyXzZ12  +  (0)Z21  +  (0)Z31  + 

(0)Z41 


zi2  “  xxXyXz  [F(F+M)(F+MKmzI>] 

2.  r(X)  -  (X+Xx)(X+Xy)(X+Xz) 

r'(X)  -  3X2  +  2(Xx+Xy+Xz)X  +  (XxXy+XxXz+XyXz) 
(F+XxI)(F+XyI)(F+XzI)  *  ^x^y^z2ll  (XxXy+XxXz+XyXz)Z^2 
+(0)Z2i  +  (0)Z31  +  (0)Z41 

Zu  -  (XxXyXz)"l|l  fJ  ( F+XxI ) ( F+Xy I ) ( F+Xz I ) 


(29) 


(30) 


3.  r(X)  -  X(X+Xy)(X+Xz) 

r’(X)  -  3X2  +  2(Xy+Xz)X  +  XyXz 

F(F+XyI)(F+XzI)  -  XyXzZi2  +  [-X3x+(Xy+Xz)X2x-XxXyXz]Z2i 

z2i  -  [-■x3x+(Xy+xz)x2x-xxXyXz]"1F Ji  -  (mxi)](F+Xyi)(mzi^i) 


9 


4.  r(\)  ■  X(X+Xx)(X+Xz) 

r'(X)  -  3X2  +  2(Xx+Xz)X+XxXz 

F(F+XxI)(F+XzI)  -  XxXzZ12+(-Xy)[X2y-(Xx+Xz)Xy+XxXz]Z31 
Z31  "  t-X3y+(Xx+Xz )X2y+XxXyXz ] "lF [( F+XxI ) ( F+Xz I )- 

(  r")(mxX  >  <  F+V  )  (  I)] 

5.  r(X)  -  X(X+Xx)(X+Xy) 

r'(X)  -  3X2  +  2(Xx+Xy)X+XxXy 

F( F+XxI ) ( F+XyI )  -  XxXyZ12  +  (-Xz)[X2z-(Xx+Xy)Xz+XxXy]Z41 
Z41  *  [-x3z+(Xx+Xy)X2z-XxXyXz]-lF[(F+XxI)(F+XyI)  - 

(fc) (F+  xI)<F+V)(mzI)]  • 

If  Equations  (29)  tnrougl.  (33)  are  now  substituted  into  Equation  (28)  the 
desired  expression  for  eFt  will  be  obtained  as 


•  Ft 


-  (x^y'x  j ( )(F+xyI)(  1  >  + 

-Xxt 

- = - 2 - - -  F(F4-XyI)(F+XzI)  + 

[-X  X+(Xy+XZ)X  X_XXXyXZ] 


-Xvt 
e  y 


[— X3y+( Xx+Xz)X2y— XxXyXz] 


F(F+XxI)(F+XzI)  + 


-X-t 


- = - = - = -  F(F+XxI)(  F+Xyl)  - 

l**x  z+(Xx+Xy)X  z“XxXyXz] 


-"Xrt 


(XxXyXz)2  ^j^y^z  l“^x+(ky+Xz)x3x-XyXzX2x] 

-Xvt 


[“X^y+( Xx+Xz )X^y— XxXzX  2y ] 


-X-t 


- T - - - = - 5—  F(  F+XxI )  (  F+Xy  I )  (  F+Xz  I ) . 

[-*V^x+VX  V^y^zU 


\yf\  z 


(34 


Equation  (34)  seems  formidable;  however,  if  specific  values  are  substituted 
for  Xx,  Xy,  Xz  the  coefficient  terms  simplify.  Equation  (34)  mainly  involves 
several  combinations  of  matrix  multiplication.  However,  for  the  general  case 
indicated  in  Equation  (34),  when  the  matrix  multiplications  are  performed  and 
the  terms  are  combined  as  Indicated,  the  resulting  exact  expression  for  the 
state  transition  matrix  (for  the  X's  of  this  section)  is  as  shown  in  Equation 


C.  Case  With  Xx,  Xy,  Xz  Positive,  Equal 

For  the  case  where  Xx  ■  Xy  ■  Xz,  the  minimal  polynomial  of  F  is  found 
to  be  1 

¥(X)  -  X2  (X+Xx)  (36) 

and  the  fundamental  formula  for  f(F),  as  given  by  Equation  (12),  is  as 
follows, 

f(F)  -  f(X1)Z]_1  +  f'(X1)Z12  +  f(X2)Z2i 

-  f(0)Zu  +  f’(0)Z12  +  f(-Xx)Z21.  (37) 

If  one  again  chooses  f(X)  *  e^fc  then  one  has  f(F)  *  eFt  so  that,  using 
Equation  (37) 

e^c  ■  Zii  +  tZ^2  +  e  X*lZ21.  (38) 

The  components  of  F  can  be  found  by  using  the  same  procedure  as  in  paragraph 
B,  with  r(X)  given  as  r(X)  »  r(Xi)Zu  +  r'(Xi)Zi2  +  r(X2)Z2^.  Proceeding 
as  in  paragraph  B  above,  one  has  the  following  results. 

1.  r(X)  -  1 
r'  (X)  -  0 

I  •  Zn  +  Z21  (39) 

2.  r(X)  -  X(X+Xx) 
r'(X)  -  2X  +  Xx 

F(F+XxI)  -  Z12  (40) 

3.  r(X)  -  X2 
r' (X)  -  2X 

F2  -  Z21  (41) 

From  steps  (1)  and  (3),  Zn  may  be  found  to  be 

ZU  -  I  -  F2.  (42) 

If  Equations  (40)  through  (42)  are  substituted  into  Equation  (38),  the 
expression  for  e“  will  be  obtained  as 
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eFt  -  I  -  F2  +  F2t  +  Ft  +  e"  x  F2 

-  I  +  (e~Xxt+t-l)F2  +  Ft  (4 

When  Equation  (43)  is  expanded,  the  resultant  matrix  expression  for  eFt  (for 
the  X's  of  this  paragraph)  is  as  follows, 

100  tOO  a  0  0 

010  OtO  0  a  0 

001  OOt  00a 


10  0  b  0  0 

0  10  0  b  0 

0  0  1  0  0  b 


c  0  0 

0  0  0  c  0 


0  0c 


where 


Wilts  Lt?  . 

a»t-l  +  ex  ( 

b  -  (1-Xx)t  +  (l-e~Xxt)Xx  ( 

c  -  a  +  (e"Xxt-l)Xx2  +  (Xx-l)Xxt.  ( 

0.  Case  With  One  Root  Zero;  Remaining  Roots  Positive,  Unequal 

To  begin  this  paragraph,  assume  that  Xy  -  0  and  that  Xx  and  X.z  are 
positive,  unequal.  Based  on  this  assumption  the  minimal  polynomial  of  F  is 
found  to  be 

?(X)  -  X3  (X+Xx)(X+Xz) .  ( 

For  the  conditions  of  this  paragraph,  the  fundamental  formula  for  f(F)  is 

f(F)  -  f(-Xx)Zu  +  f(0)Z2i  +  f  ’  (0)Z22  +  f"(0)Z23  +  f(-Xz)Z31  ( 

and  eFt  will  be  given  by 

eF^  m  e  Xx^Zn  +  Z2i  +  tZ22  +  t2Z23  +  e  Xz^Z3^.  ( 

The  components  of  F  are  again  found  by  the  same  procedure  used  in 
paragraphs  B  and  C  where,  in  this  paragraph,  r(X)  is  given  by 
r(X)  -  r(-Xx)ZU  +  r(0)Z21  +  r’(0)Z22  +  r"(0)Z23  +  r(-Xz)Z31. 

1.  r(X)  -  X3(X+Xx) 


1 


ituted  into  Equ. 
paragraph)  Is  g 


10  0 
0  10 
0  0  1 


-aX  3x+bX  2  t^+t 


aX2z-bXz 


0 


-aX3z+bX2z-l/2Xzt2+t 


X-e  ^xt  /  X3,+X2-Xr+XvX2,+X3~\  / X^+X-yXy+X2 


b  *  -  .  + 

^  x(^x"^z) 


/  Xx+^z  V  2  _ 

V  2\X\J 


^  x  ^  x/vz^'vx/v  z  ^  z 
(  X3xX3z 

i  ~xzt 
Xxe  z 


x^x^z™  z 

^V2Z  ) 


^■^z(^x~^z) 


For  a  case  where  Xx  -  0,  with  Xy,  Xz  positive,  unequal,  the  minimal 
polynomial  of  F  is 

¥(F)-  X3(X+Xy)(X+Xz)  (57 

and  the  fundamental  formula  for  f(F)  is  given  by 

f(F)  -  f(0)ZU  +  f’(0)Z12  +f"(0)Zi3  +  f(-Xy)Z2l  +  f(~Xz)Z31  .  (58 

In  a  similar  fashion,  for  Xz  -  0  and  Xx,  Xy  positive,  unequal,  the  minimal 
polynomial  of  F  is 

¥(X)  -  X3(X+Xx)(X+Xy)  (59 

and  the  fundamental  formula  for  f(F)  is  given  by 

f(P)  "  f(“^x)Zll  +  f(-^y)Z21  +  f(0)Z31  +  f * (0)Z32  +  f"(0)Z33  .  (60 

The  expressions  for  e^c  corresponding  to  Equations  (58)  and  (60)  are  given  in 
Equations  (61)  and  (62),  respectively. 


L5 


2 


X2x+XxXy+X2. 


X2x+XxXv+X2 


E.  Case  With  One  Root  Zero;  Remaining  Roots  Positive,  Equal 

For  a  case  with,  for  instance,  Xx  -  0,  Xy  ■  X2  >  0,  the  minimal  poly 
nomial  of  F  is  found  to  be 

¥(X)  -  X3(X+Xy)  (63 

and  the  fundamental  formula  for  f(F)  is  thus  given  by 

f(F)  -  f(0)ZU  +  f’(0)Zi2  +  f"(0)Z13  +  f(-Xy)Z2i.  (64 

The  expression  for  eFt  is,  therefore,  given  as 

eFt  m  Zii  +  tZ^2  +  t^Zi3  +  e  ^^21  (65 

and  the  components  of  F  can  be  found  as  follows: 

1.  r(X)  -  X3 

r'(X)  -  3X2  ;  r"(X)  -  6X 

F3  -  -X3yZ2i 

Zn  ' "  f3  <m 

2.  r(X)  -  X2(X+Xy) 

r'(X)  -  3X2  +  2XyX  ;  r"(X)  -  6X  +  2Xy 

F2(F+XyI)  -  2XyZ13 

zi3  ■  (s;)f2<p*V>  <67 

3.  r(X)  -  1 


r* (X)  -  r"(X) 
I  -  Zu  +  Z2i 
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IV.  CONCLUSIONS 


In  order  to  verify  the  expressions  developed  for  eFt  in  Section  III,  a 
computer  program  was  written  which  would  integrate  the  target  equations  (2), 
for  an  arbitrary  set  of  initial  conditions  and  X.'s  and  for  w  =  0,  to  generate 
time  histories  of  target  position,  velocity  and  acceleration  along  the  x,  y, 
z  axes.  A  listing  of  the  program  is  presented  in  Appendix  A.  Using  this 
program,  data  was  generated  for  each  of  the  cases  considered  in  Section  III. 
Table  1  lists  the  target  state  initial  conditions  and  the  values  used  for  \x , 
Xy  and  \z  for  each  computer  run  case.  The  target  state  transition  matrix  from 
Section  III,  which  corresponds  to  each  of  the  cases  listed  in  Table  1,  was 
then  used  to  generate  target  state  data  for  comparison  with  the  computer  out¬ 
put. 

A  representative  set  of  data  indicating  the  results  of  the  comparison  is 
presented  in  Figures  1  through  12.  This  data  represents  the  target  position, 
for  each  of  the  cases  in  Table  1,  for  both  the  computer  solution  and  the  state 
transition  matrix  solution.  As  can  be  seen,  the  data  output  matches  in  each 
instance  and  the  state  transition  matrix  expressions  in  Section  III  are 
verified . 


TABLE  1.  TEST  RUN  CONDITIONS 


Paragraph 


NOTE:  Target  state  initial  conditions  in  each  case  were 
as  follows: 


STATE  TRANSITION  MATRIX 


Figure  1.  Target  x  Position  Versus  Time,  Case 
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Figure  4.  Target  x  Position  Versus  Time,  Case 


»SE  3:  TIME  VS.  POSITION 

ATE  TRANSITION  MATRIX  -  A  29  MOV-02 


Position  Versus 


5:  TIME  VS.  POSITION 

TRANSITION  MATRIX  -  A  29- WO V- 02 

ID :  24: 43 


Figure  11.  Target  y  Position  Versus  Time,  Case 


Figure  12.  Target  z  Position  Versus  Time,  Case 


APPENDIX  A 


DIGITAL  SIMULATION  PROGRAM 


A  listing  of  the  digital  simulation  program  which  was  used  to  integrate 
the  differential  equations  describing  the  target  acceleration  components. 
Equation  (2),  in  order  to  provide  time  histories  of  target  position,  velocity 
and  acceleration  components  is  presented  in  this  Appendix.  The  simulation 
permits  target  state  initial  conditions  and  values  for  \x,  Xy ,  \z  to  be  varied 
by  the  user  and  uses  an  external  line-printer  plot  routine  to  provide  data 
plot  outputs. 


-  L - 1  i  *1-  PROGRA/*  WiLl  T 

;v  ■-  —  1  .  ‘  aP\tL  :  Av_  Cl  z.H  AT  I  tjN'E  ‘.El!1*’.' 

i-‘£  FR  .'Gw.-M  will  USE  ~hE  RuM&E-i-.UT- 
.  1  4  •  .  C  x  •  9  • .  # UT Ta.  DT  .  N  >. 

^  -  .  .  .5  ; d/5.  s  a.  i . 

•-••-L.  •••  -R.  IW,  15.  ID; 

v  .;.l  *4  '  - 

'  Ox  n  r  i  -  . 

OU‘  ;  *  -  10 

I  I 

■_  f  i  T  i  :-ii. 

*  -  0 

RfeAD  i “4 1  ;  ai_  position,  velqc i  r  <  ac ce 

~'Fau  »  IP.  «•  t£ND=24 1  <  <  3  » •  X  « 2  ) ,  a  ;  x  > 

aE  AO  ■  IP-  «■  i  <  •.  6  )  ■  X  <  5  !  .  x  ■  4  ) 

At  AD  .  I P  ■  *  '  X  <  9  >  •  X  \  S )  •  <  7  j 
~  r  E  I- 

''j  /:  =  4 

READ  uAMSBa  max.  ay  xkZ 

READ 1  IP-*  »n»  XKV,  X*Z  ; 

-ONTINoE 

IF  s  T I fi£  •;£  20  0 )  GO  TO  2 1 

*  *  ft  ♦  1 

DO  1 OC  *v  T  f  a  =  i,4 

D  X  '  1  (  =  <  -  *  A  A  <  fr  X  .  ;  > 

0  i  i  2  i  =  <  ■:• 

D  X  ( 3  j  -  »  2 » 

DX  i  4  ; 

D  X  t  5  i 
D  X  i  6  . 

jX-7;  =  l-  ««.2i*X,7) 

D  x  <  9 1  *  x  t  > 

D  X  <  9  >  *  X  ■■  S  • 

GO  TO  <  3o. *0. 00.  40/ . PUTT* 


=:  t  —  .•  a  v  i  *  x  '  4  J 
=  X  i  4  ■ 


CONTINUE 

IF  <K0Ty>  GT  i  yu  Tv  01 
CONTINUE 

T  iME  =  Ti  '»E  -  v  *CT 

CONTINUE 

CALL  ROnv 

IF  ( ftUUN'i  r.E  0-  GO  TO  C4 

WRITE  ■  6-  214)  TIME,  X  <  i  .  u  2  .»  •  x  j; 

FORMAT!  2X-  4..012  o- 2X  . 

WRITE'  6 . 2 1 o . X  >  4  . ,  x ■ 5 .  «  6> 

wRIT£ It,.  ria. ,  •  -  x  '3 
FQRma  r  •  2  "■  J  •  G  i  2  o  ■  2  X 
C  DN~  I  i'l'.-'E 

*•>  M  S  »♦  ■*  ’ 

■  -  •  *  0  jN  "  £  j  *r  0  v  W T 


=  C  "-'AEt  Cj'jA 

I  TIES.  A-;£  P!,;  f  t-  .OMO 
A  4  SuBhv'  r INE 


l£«AT  iuf.  FOR  x  v 


/ 


.  f  i : 


WRITE  (to, 200 i TIME,  X(1>,X<2).X<3> 
FORMAT <2X,  4<G12.  6,  2X> ) 

WRITE  <6, 205)Xl4),  X<5),  X<6) 

WRITE <6, 205) X< 7), X<8>. X<9) 

FORMAT  (2X, 3(G12. 6, 2X) ) 

CONTINUE 

KOUNT 1 =KOUNT 1+1 

IF (KOUNT 1  EQ.  MKOUNT )  R0UNT1=0 

GO  TO  1000 

CONTINUE 

GO  TO  511 

STOP 

END 


SUBROUTINE  RUNK 

COMMON  X  <  9) , DX  <  9 ) , KUTTA,  DT ,  NX 

DIMENSION  XA<9), DXA<9) 

GO  TO  (10,30,50.70),  KUTTA 
DO  20  I  =  1 ,  NX 
XA< I )  =  X ( I  ) 

DXA < I )  =  DT#  DX( I ) 

X  <  I  )  »  X  (  I  )  +  5*DXA  <  I  ) 

RETURN 

TDT  =2  *  DT 

HDT  *  .  5  #  DT 
DO  40  I  =  1 , NX 

DXA(I)  -  DXA < I )  +  TDT  *  DX < I ) 

X  C I  )  -  XA<  I )  «■  HDT  *  DX  <  I ) 

RETURN 

DO  60  I  *  1 , NX 
VDT  =  DT  *  DX ( I ) 

DXA ( I )  =  DXA(I)  +  2.  *  VDT 

X ( I )  *  XA< I )  +  VDT 
RETURN 

DO  80  I  =  1 ,  NX 

X(I)  =  XA  <  I  )  +  (DXA(I)  ♦  DT  *  OX  (I,-  i/6 
RETURN 
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